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Proc6d* de determination automatique du contraste et de la 
brillance d'une image radiographique numerlque 



^invention concernc la determination automatique de la 
bnllance et du con.raste dune image radiographique numeric d'un 

0t>jet ' EUe trouve une application particulierement interessante dans 
,c domaine medical dans lequel on proc.de a des analyses fines de. 
in.agesradiogr.phiques.enpaniculierlamanunograph.e. 

Dans le domaine medical, les diagnoses se basent 
eendralement sur l'etude dune image radiograph.de. Le progres des 
"um^nques fait q.il -I aujourdh. possible de faire v^r 
les caractenstiques do Vimage afm de realiser au nueu* le diagno^. 
U est done evident que la qualite de 1'image est un point essenuel. 

Comme cela est bien connu. lacquisition d image fait 
intervenir de nombreux paramctrcs Hes a la chaine ^acquisition : on y 
di "gue les paramo ^acquisition (matenau ciblc de 1 anode du 
be rayons X (piste), materiau et epaisseur du nitre en sortie d 
c a rayons X (filuc). -ion appliquee entre ca^ ^ j 
tub e a rayons X (kV). courant de chauffage du filament de ,a cathode 
do tube a rayon. X (mA). temps deposition. ...). les P™ 
positionneur (incidence de la prise de vue, facteur dagmndt seme,, 
tvpe de pelote de compression utilisee. epaisseur de lobjet for, de 
compresL, ) et les parametre, du detecteur digit* ,el auon 
niveaux de gris-flux capte de rayons X....). F.n plus de ce.a un Lu 
fctervenir les parang de 1'objet. notamment sa composition. Ca 
ensemble de paramfctres P ermct d'obtenir une image qui est en lai 
compos* de pixels plus oo moins grises. La quantification de . image 
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est obtenue sur un nombre donn£ de niveaux de gris. Lorsquc l'imagc 
est affichee sur un ecran ou imprime sur un film en utilisant unc 
gamme de niveaux de gris donne\ le contrastc percu peut 6tre 
considdrablement different entre les images du mfime objet acquiscs 
avec differents ensembles de parametres d'acquisition ou entre des 
images de differents objets obtenues avec le mime ensemble de 
parametres. En outre, le choix de la gamme de niveau de gris a 
considercr pour l'ecran est d'une grande importance afm de donner un 
percu de contraste acceptable dans l'image. Ainsi, une fois l'image 
obtenue, l'utilisatcur modifie interactivement la brillance (WL : 
window level) et le contraste (WW : window width) de cctte image 
afin d'y ajuster la gamme dc niveaux de gris de facon a atteindre un 
percu de contraste correct. 

Generalement pour plus de 16000 niveaux de gris different*, 
il est difficile et long dc choisir manuellemcnt la brillance ct le 

contrastc corrects. 

L'invention vise a apporter une solution a ce problemc. 

Un but de l'invention est de determiner automatiquement la 
valeur de la brillance (WL) ct d'en deduire le contraste (WW). 

L'invention propose done un proc6de d'acquisition d'imagc 
radiographique numenque d'un objet avec un reglage automatique des 
parametres de visualisation (brillance, contraste). Ce precede est 
adaptable, puisque les scuhaits de l'utilisatcur peuvent ette pris en 
compte. 

Ce proceed comprend une obtention de deux niveaux de gris, 
gris min et gris max revelant unc zone dc niveaux dc gris dans laquelle 
la brillance (WL) est susceptible d'etre determined un jeu de regies 
permet le calcul dc cette valeur WL. et on en deduit fmalemcnt une 
gamme dc niveaux de gris pour le. contraste WW. 

Selon une caracteristique geneVale de l'invention. on effectue 
une premiere elape dc segmentation. Les donnees d'entree de cctte 
premiere etapc sont les parametres d'acquisition. le. paramet.es du 
ddtecteur digital, (es parametres du positionnemcnt et les parametres 
de l'objet L'cnsemble de ces parametres est introdmt d,ins un modele 
mathCmattque prf-etabli de la chainc image et de l'objet. 
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Ce modele permet de determiner un niveau dc gris mm et un 
niveau de gris max. La partie encadr6e par ces deux niveaux est la 
zone utile des niveaux de gris dans laquelle il sera choisi la valeur du 
niveau de brillance WL, 

Cette £tape correspond a une segmentation de l'image par un 
double seuillage de telle sorte que soils les points de l'image dont le 
niveau de gris min et le niveau de gris max soient conserves. 

Une seconde dtape est realist parallelcment a la precedence.' 
Elle met en oeuvre la veritable acquisition radiograph ique de I'objet. 
A partir de l'image obtenue, un histogramme est dresse : l'histogramme 
r<Sel. 

La troisieme etape fait intervenir gris min. gris max et 
l'histogramme reel. Leur mise en correspondance abouiit a la 
suppression, dans l'histogramme reel de la partie inferieure au niveau 
de gris min ct de la partie superieure au niveau de gris max. Le 
resultat de la troisieme etape est un intervalle de niveaux de gris 
limitant la zone utile : histogramme limile. II reste done a determiner 
la valeur de la brillance WL a I'intdrieur de cet intervalle. 

La quatrieme etape conceme done la determination de WL 
qui peut se faire au moyen dun jeu de regies simples connu de 
l'homme du metier Ce jcu de rfegles permet d'obtenir une valcur WL. 
Cependam lorsque rutilisateur precede a de nombreuscs 
mammographies, il est n6cessaire d'ameliorer le rendement en faisant 
intervenir un modele mathematique de l'histogramme. Un modele est 
decrit par M. Jean Lienard LSD/AAP reference 98030 Technical Note 
GE Medical Systems. August 1998. II considere un sein comprint 
dont la forme est modelisee, par example, a l'aide de formes 
geometriques simples. En faisant intervenir des donnas du sein a 
mammographier, on obtient un modele d'histogramme caractense par 
un maximum au niveau de 1'epaisseur radiologique corrcspondant a un 
tissu du sein. 

La mise en correspondance du modele mathematique do 
rhistogramme ainsi determine el de l'histogramme hmite permet de 
determiner la valeur d'une epaisseur radiologique caractdrisant le sein. 
35 cet histogramme Hmite ainsi affine va subir un jeu dc regie. 
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permettant la determination dc la brillance WL. 

La cinquieme etape .oncerne la determination du contraste 
WW Cette valeur pent etre obtenue d'une manure connue par l'hommc 
du metier en faisant mtervenir la valeur WL et d eventuels paramfetres 
connus a prion ct dependants dc 1'utilisatenr. 

Selon un mode de mise en oeuvre de Invention, le contraste 
WW pent egalement etre determine ind^pendamment de WL. On utilise 
dans ce cas le coefficient d'attenuation moyen (jxmoyen) dependant du 
spectre. D'une facon gentle, le WW ainsi obtenu est fonct.cn du 
spectre, car proportionnel a \i moyen. 

LoLn.ion de WL. et WW (.ndependamment de WL) eqmvau, 
en fait a une operation d'au.o-con.raste pouvan. e.re util.se dans un 
precede de perception des centres en epaisseur invanante. 

Un but de ce procede e S . de rendre la percept.on d es 
contrastes independante des condi.ions ^acquisition e. de , r«*j«. _Pou. 
deux objets d'epaisseurs different on veut que que.que ou le 
conditJs d'acquisi.ion, le contraste percu sur . itna e e teste f.dele 
lew difference d'epaisseurs recllcs. ^ 
En d'autres .ermes. on vcut qu'une epa.sseur donne reprfse, 
roujours un meme per !U de contras.e quelque soi. lobje, b. 

conditions dacquisition. /Uments 
En effet au cours d'une acquis.tion a. mage, les elements 
piques fon, que ,e specue transform 1 epaisseur reelle de roojet en 
p ,S eur rad,olo g ique entrainan, a,nsi un percu d= con.ras e .ntuWe. 
Ces. un effe. dat.enuation exponent.clle. On compense cc, e.ie. en 
■nLuisant sur une c haS ne ima S c. un etage d= clement espae, 
afl „ dealer Vat.cnua.ion exponent.elle due a l .n.eracon des rayon 
X avee lobjet cn employam une function logan.hmique mod.f.ee. Le 
cUngemen, d'espace pe.net de quiuer fespacc exponent,, pot,, 
passer dans lespace des cpa,sseurs rad.olog.ques. 

On imrodui. cgulement un 6ta e c de visuahsauo,, pcnnet.an. 
de passer de Vespace des cpaisseurs rad.ologiques vers lespa« de, 
eparsseurs reelles. Ce chan^n, despacc «. poss.blc .a U W e 
proportionnel a U coefficient d'at.enuauon Uneaire moyen de .obj... 
etudi£ ct independant de WL. 
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D'autrcs avantages et caracterisUques de Hnvent.on 
apparattront a Vexamen de la description detaillce d un mode de mtsc 
en oeuvre nallement linritatif, et des dessins annexes sur lesqueh : 

- la Figure 1 est un orgaiugramme d'un mode de mise oeuvre 

S du procede selon 1'invention, 

- la figure 2 illustre l'histogramme limit* par rapport a 

l'histogramme reel, 

- la figure 3 illustre une fa^on de choisir WL, 

- la figure 4 decrit le procede detention d'un coniraste 

10 ^"^Bien que Invention ne soit pa, Umite,, on va nontenant 
decnre location du precede salon invention a ^munat^ 
automate de la brillance WI, et du contrastc W* pour 

^n.ere -pe met en oeuvre un mod.e mad.m^ue de 
la chaine image et de 1'objet avec en donnee. d'entree les parametres 
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suivants : 



'. epaisseur d- sein comprint et parametres du positionneur 
(incidence de ,a prise de vue. fac.eur dagrandissement. type de peloK 
de expression utiUsee, epa,sseur de r«**«. force de -P— ■ 
..)corrmteparametresdu positionneur, 

. parame.res du detecteur (relation en.re le flux de rayons X 
recu sur le deteeteur e. les niveaux de gris de i rmage produ.te . . ) 
- parametres d'acquisition (piste, fillre. kV. mAs, ...) 
. papametres d'obje. (epaisseur mccanique du S e.n. valeu.s 
minimum um.net maximum Umax du coefficient donation hnauit 

de 1'objet, ^ ^ pri^ijalOTem compose de tissus fibre* et 
,dineux si Ton na aucune information sur la con.porit.on du se.n, on 
enitaer ,m,n e, .max en faisan, deu< suppos.tions extreme, 

p.min correspond „ coefficient dattenuation ,ne.,.r dc 
ti , sus les moms attendants de lobje, (.i»»s ad.peux pour le sent) pou 
ies energies du spectre de rayons X determine par les pa.amet.es 
d'acquisition. 
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Umax peut 6tre estime" de deux manieres : d'une part commc 
u.min, en eansiderani que Umax correspond au coefficient d'attenuanon 
lineaire des tissus les plus absorbants de 1'objet (fibres pour le sein) 
pour les Energies du spectre de rayons X determine par les parametres 
d'acquisition ; d'autre part d'une maniere plus precise a partir dun 
modele mathematique de la chaine image, de 1'cpaisseur mecamque du 
sein comprint des parametres d'acquisitlon et dune quanttte de 
photons obtenue a la suite d'une propose realisee sur une zone de 
Lsite maxtmum (qui permet d'est.mer la valeur du coetf.cten 
d-attenuation lineaire corresponds a la zone la plus attcnuante de 

r ° bjCt> " L'ensemble de parametres iniroduit dans le model* 
m athemati q ue de la chaine image et de ,'objct permet d'avoir a .a 
sortie deux valeurs de niveaux de gris : gris min et gns max (hg. 
Ces deux valeurs delimiter* en fait la zone utile qui est unc zone de 
niveaux de gris se rapportant reellement au sein. En ™ 
deux valeurs extremes de la composition du sem ^ ct on 
obtenu une zone delta**. par deux valeurs e.^mes gns m n et gr 
ma x en dehors de Quelle les niveaux de g ris ne pa> - 

sen. Plus preelaemenu la panic de niveaux de gm.nf ncurs - a r 
min correspond a des objets plus attenuants que robjet ~ 
partte de niveaux de gns superiors a gris max correspond au fond d 
Lege. Cette etapc est unc etape de ^mentation car elle permet de 

delimiter la zone utile* j^ 1v 

On cflfecuc aiors une mi s e en «■« '« 

vaieur, gns min e, gris - « "» hi s ,o g ramm= tau * I .muge 
;i„ g rap hiqU e du se.n (nts,o g « re>l, Plus pre— . . 
«W. to panic injure a gris m,n e, la P-« «•«— = 
m „ de Jon a nc cc.nser.er <,,e la ™ne u.«c : histograms l.m.te 

( '"' g La brdlance WL «. une v.ieur toctee to »»« «• 

pent etre obtenue de plusieurs facons. 

Une maniere dobtentton de WL est lapphcauon d un jcu 
rig ,cs prt-etablies a Ito.ogramme limi... Un jcu de regies P«" 

comprendre : 
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- determination du niveau de gris correspondant au maximum 

dc 1'histogramme limits 

- conservation d'une quantite de \% (typiqucmcnt 95tt) des 
occurrences dudit histogrammc limit* a droite du maximum, ct x% 
egalemcnt de S occurrences a gauche du maximum : on obtient ainsi un 
histogramme reconstruit 

- la determination de WL en tam que valeur mediane de 

1'histogramme reconstruit (fig- 3)- 

U C st egalement possible d'obtenir WL avec une me.lleure 
precision en faisant preceder 1'etape de jeu de regies dune etape de 
mise en correspondance. Cette etape fait intervenir un modele 
mathematique de rhistogramme dans lequel : 

- la forme du sein est un cyl.ndre genere par rotation autour 
dun axe, d'un rectangle dent Vun de, petits c6tes est ferme par un 
demi-cerele de diametre 6gal a la longueur de ce pent c6te. 

- la composition du sein est homogene, par exemple 100 de 

EiaiSSe - on etablit un histogramme qui correspond a la probability 

en fonction de l'epaisseur radiologique du sein r<(|Vli « cu , 

- le maximum de rhistogramme obtenu represente lepaisscm 
rn~ du sein qui, multiple par lc coefficient dattenuation 
donne 1'epaisseur radiologique maximum correspondant au tissu 

adipeux (graisse). , 
La mise en correspondance des deux histogrammes (modele 

mathematique « histogram** limit© pcrmet de ^—J^ 
1-histogranune limite la valeur de lepaisseur radiologique maxnnum 
correspondant au tissu adipeux dans le sein. 

Ccci pcrmet done de determiner les valeurs dcpa.sseur 
radiologique des differents composants du sein. 

Cette mise en correspondance est reaHsee en employant une 
methode do minimisation d'errcurs euUc deux fonctions comme par 
cxcrnple la methode des rnoindres carres. 

On peut ensuite appliquer un jeu dc regies adequat afin de 
determiner WL. Par exemple WL=aE avec E represcntant la valeur dc 
1-epaisscur radiologique obtenue a partir du modele mathem.t.qu, dc 
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l'histograirirne. 

L'6tape suivantc correspond a la determination de WW en 
utilisant WL. On obtient WW a partir d'une fonction faisant intcrven.r 
done WL et eventuellement d'autres parametres, notamment un 
5 parametre G (fig. 3) qui est choisi par l'utilisateur. 

Ce parametre G rend done ce precede adaptable a chaque 

utilisateur 

WW = g(WL, G) 

g etant une fonction qui a partir de WL et G determine WW 
10 qui represente done une gamme de niveau de gris autour de WL 

II est egalement possible de determiner WW independamment 

dc WL. . , 

On determine d'abord u.moyen a partir des informations tnees 

de I'histogramme limit*. Comme exemple, on pent prendre le M 
15 correspondant a la vleur mediane de I'histogramme limit* comme 

umoyetv _ 
Ensuite, on utilise une loi faisant intervenir une constants 

Cte afin d'en d6duire WW : 

WW = Cte \i (spectre) 

Cette relation est vra.e dans un cas monoenergehque mats 
d'une facon generate, WW est fonction du spectre 
WW = f (spectre) 

L'operation d'auto contraste est ainsi realisee puisquc WL ct 

WW sont etablis. 

On va niaimenant decrire la facon dont la percept.on du 
central d une difference d'epaisscur reste invariante quels que scent 
les moyens d'acquisition et Its objets conformement a la figure 4 

^intervention des rayons X avec lobjet cntrainc unc 
attenuation exponentielle de lintensite 1. an niveau de l'image : 
%Q 1 = lo exp (-Wl) 

lo e<t une constant* 1 represente une grandeur infim.esimale 
qui correspond a une distance le long du trajet reliant le foyer da 

rayon X ct le detecteur. 

1 udl represente l'epaisseur radiologique pour une z 

35 donn6e d un objet. 
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C'est cette quantite qui nous int£resse. Pour l'obtenir, on 
realise une operation de prf-lut (look-up table) permettant de 
compenser l'attenuation exponentielle en utilisant une fonction 
logarithinique modified. Elle est dite modifiee car les niveaux de gris 
S les plus bas sont transformes scion une rotation lindaire alors que la 
fonction logarithmique intervient progressivement pour les autre s 
niveaux de gris. 

On se trouve ainsi dans l'espace des 6paisseurs radiologiques 
dans lequel une epaisseur radiologique peut etre notee Jill avec H 

10 l'epaisseur reelle. 

On veut finalemcnt realiser un changement d'espacc pour sc 

retrouver dans l'espace des dpaisseurs belles. 

L'autocontrastc nous foumit les valeurs WW et WL que 1'on 
introduit au niveau d'une lut de visualisation. Ladite lut de 
visualisation pcrmet d'eliminer u.. Dans le cas dune image 
monodnergetique. cette operation revient a une division par u. 

Cette operation est possible car WW est proportionnel a u. 
Le nSsultat ainsi obtenu pcut etre introduit dans le systeme de 
visual.sation standard DICOM (Grayscale Standard Display Function, 
suppldment 28) connu par rhomme du metier afin de visualiser Hmage 
(fig- 4). 
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REVENDICATIONS 



1. Precede d'acquisition d'une image radiographique 
numenque d'un objct. comprenant un histogramme de Vimage 
radiographi6e (histogramme reel), un modele mathematique de la chat- 
ne image et de l'objet obtenu par calibration, caracterise par le fait 
que: 

a) on se sert du modele mathematique de la chaine image ct 
de l'objet et dun ensemble de parametres d'acquisition, du detecteur, 
du posttionneur el d'objet pour determiner deux valeurs de niveau de 
gris, gris min et gris max delimitant une zone utile, 

b) on supprime dans l'histogramme reel la partie infdricure a 
gris min et la partie superieure a gris max : histogramme limitd, 

c) on applique un jeu de regies audit histogramme limite afm 
de determiner le niveau de brillancc WL, 

d) on obtient le contraste WW soil a partir de WL ct 
eventuellement d'un ou plusieurs paramfctrcs choisis par I'utilisateur ou 
fix6s a priori, soil a partir d'un coefficient d'attenuation moyen de 
l'image indfpendammeni de WL. 

2. Procdde selon la revendication 1, caraciensC par le fait que 
le coefficient d'attenuation lineaire minimum umin de l'objet est 
estime a partir de valeurs connues du coefficient d'attenuation lindane 
des tissus les moins attenuants de l'objet (tissus adipeux pour le sien) 
pour les energies du spectre de rayons X determine par les parametres 
d'acquisition et permct avec les parametres du detecteur, les 
parametres du posiUonneur, les parametres d'acquisition et les 
parametres de l'objet, au travers d'un modele mathematique de la 
chaSne image ct de lobjet, dc determiner gris max. 

3. Procede selon la revendication 1. caractense par le fait que 
le coefficient d'attenuation lineaire maximum (umax) de lobjet est 
estime soit a partir dc valeurs connues du coefficient d'attenuation 
lineaire des tissus les plus au6nuants de l'objet (tissus fibreux pour le 
sein) pour les energies du spectre dc rayons X diamine par les 
parametres d'acquisition, soit de manierc plus precise a partir dc 
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l'epaisseur mecanique du sein comprint, des parametres d'acquisition, 
du modele mathdmaliquc de la chame image et de l'objet et d'une 
quantity de photons obtenue h la suite d'une propose realised sur une 
zone de densite maximum de l'objet, et pcrmet avec les parametres du 
detecteur, les parametres du positionneur, les parametres d'acquisition 
et les parametres de l'objet, au travers d'un modele mathematique de la 
chalne image et de l'objet, de determiner gris min. 

4. Proc^de scion l'une quelconque des revendications 
precedentes, caracterise' par le fait que la limitation de l'histogramme 
fait intervenir deux niveaux de gris (gris min et gris max) encadrant 
une zone de niveaux de gris dans laquclle la brillance WL est 
d6termin.ee. 

5. Precede" selon l'unc quelconque des revendications 
precedentes, caracterise par le fait qu'on realise une mise en 
correspondance cntre l'histogramme limits et un modele mathematique 
de l'histogramme afin de determiner la valeur d'une epaisseur 
radiologique caracterisant l'objet. 

6. ProcedS selon la revendication 5, caracterise par le fait que 
la mise en correspondance entre rhistogramme limits et le modele 
mathematique de lliisiogramme s'effectue en applicant une methode 
de minimisation d'erreurs entre deux fonctions. 

7. Proc6d6 selon l'une quelconque des revendications 
precedentes, caracterise par le fait que lorsque le contraste WW est 
obtenu independamment de la brillance WL, on introduit les valeurs 
WW et WL dans une chame d'acquisition comportant un eiage (Pre- 
lut) de compensation de l'effet d'attdnuation exponentielle du 
rayonnement en utilisant une fonction logarithmique modififie, de sorte 
que la perception du contraste dune difference d'epaisseur donnee 
restc invariante quclques soient les moyens d'acquisition. 

8. Precede scion la revendication 7, caracterise par le fait que 
WL et WW interviennent sur un etagc de visualisation de telle sorte 
que les Signaux sortams dudit etage (Pre-lut) de compensation (espace 
des epaisseurs radiologiques) sunt dimensionnes dans I'espace des 
epaisseurs replies. 
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FIG.1 




LIMITATION DE L'HiSTOGRAMME 
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FIG.2 



OCCURRENCES 




FIG.3 



OCCURRENCES 



WW 




WL VALEUR MEDIANE DE 
L'HISTOGRAMME 
RECONSTRUIT 

WW : g (WL, G) 



NIVEAUX DE ORIS 
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